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By treating the suicide as a social fact, Durkheim envisaged that suicide rates should be determined by the 
connections between people and society. Under the same framework, he considered that crime is bound up 
with the fundamental conditions of all social life. The social effect on the occurrence of homicides has been 
previously substantiated, and confirmed here, in terms of a superlinear scaling relation: by doubling the 
population of a Brazilian city results in an average increment of 135% in the number of homicides, rather 
than the expected isometric increase of 100%, as found, for example, for the mortality due to car crashes. 
Here we present statistical signs of the social influence on the suicide occurrence in cities. Differently from 
homicides (superlinear) and fatal events in car crashes (isometric), we find sublinear scaling behavior 
between the number of suicides and city population, with allometric power-law exponents, p = 0.84 ± 0.02 
and 0.87 ± 0.01, for all cities in Brazil and US counties, respectively. Also for suicides in US, but using the 
Metropolitan Statistical Areas (MSAs), we obtain P = 0.88 ±0.01. 

It is not uncommon in nature to observe properties that present non-trivial forms of scale dependence. This is 
the case, for instance, of critical phenomena, where scaling invariance, universal properties and renormaliza- 
tion concepts constitute the theoretical framework of a well-established field in physics^. In biology, the so- 
called allometric relations certainly represent outstanding examples of natural scaling laws. Precisely, allometry 
implies the use of power-laws, Y oc M^, to describe the dependence of a wide range of anatomical, physiological 
and behavioral properties, denoted here as 7, on the size or the body mass of different animal species, M. If the 
scaling exponent is p = 1, the variables 7 and M are trivially proportional, and the relation between them is said to 
be isometric, while j5 1 indicates an allometric type of relationship. The so-called "three-quarters law" or 
Kleiber's law, as originally proposed by the agricultural biologist Max Kleiber in 1947^, is surely one of the most 
prominent allometric relations found in natural sciences. Based on an extensive set of experimental data, this 
fundamental law states that the metabolic rate of all animals should scale to the 3/4 power of their corresponding 
masses^'^. 

In analogy with biological scaling laws, Bettencourt et aV showed that, regardless the enormous complexity 
and diversity of human behavior and the extraordinary geographic variability of urban settlements, cities belong- 
ing to the same urban system obey pervasive allometric relations with population size, therefore exhibiting 
nonextensive rates of innovation, wealth creation, patterns of consumption, human social behavior, and several 
other properties related to the urban infrastructure. The authors then conclude that all data can be grouped into 
three categories, namely, material infrastructure, individual human needs, and patterns of social activity. Despite 
the unambiguous presence of power-laws, the urban indicators do not necessarily follow a universal behavior. 
Instead the results can be divided in three different classes. The isometric (linear) case {{i=l) typically reflects the 
scaling of individual human needs, like the number of jobs, houses, and water consumption. As in biology, the 
allometric sublinear behavior (j5 < 1) implies an economy of scale in the quantity of interest, because lis per capita 
measurement decreases with population size. In the case of cities, this is materialized, for example, in the number 
of gasoline stations, the total length of electrical cables, and the road surface (material and infrastructure). The 
case of superlinear allometry (jS > 1) in urban indicators emerges whenever the complex patterns of social activity 
have significant influence. Wages, income, growth domestic product, bank deposits, as well as rates of invention, 
measured by new patents and employment in creative sectors, all display a superlinear increase with population 
size^. While these results indicate that larger cities are associated with higher levels of human productivity and 
quality of life, superlinear scaling can also characterize negative urban scenarios, such as the prospect of living 
costs, crime rates, pollution and disease incidence^'^'^"^^. 
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Figure 1 | Scaling relations for homicides, traffic accidents, and suicides for the year of 2009 in Brazil. The small circles show the total number of deaths 
by (a) homicides (red), (b) traffic accidents (blue), and (c) suicides (green) vs the population of each city. Each graph represents only one urban indicator, 
and the solid gray line indicate the best fit for a power-law relation, using OLS regression, between the average total number of deaths and the city size 
(population) . To reduce the fluctuations we also performed a Nadaraya-Watson kernel regression^^'^^ The dashed lines show the 95% confidence band for 
the Nadaraya-Watson kernel regression. The ordinary least-squares (OLS)^*^ fit to the Nadaraya-Watson kernel regression applied to the data on 
homicides in (a) reveals an allometric exponent P = 1.24 ± 0.01, with a 95% confidence interval estimated by bootstrap. This is compatible with previous 
results obtained for U.S.^ that also indicate a super-linear scaling relation with population and an exponent P = 1.16. Using the same procedure, we find 
P = 0.99 ± 0.02 and 0.84 ± 0.02 for the numbers of deaths in traffic accidents (b) and suicides (c), respectively. The values of the Pearson correlation 
coefficients p associated with these scaling relations are shown in each plot. This non-linear behavior observed for homicides and suicides certainly reflects 
the complexity of human social relations and strongly suggests that the the topology of the social network plays an important role on the rate of these 
events, (d) The solid lines show the Nadaraya-Watson kernel regression rate of deaths (total number of deaths divided by the population of a city) for each 
urban indicator, namely, homicides (red), traffic accidents (blue), and suicides (green). The dashed lines represent the 95% confidence bands. While the 
rate of fatal traffic accidents remains approximately invariant, the rate of homicides systematically increases, and the rate of suicides decreases with 
population. 



Methods 

Brazilian Data. We analyzed data available for all Brazilian cities from 1992 to 2009, 
made freely available by the Brazil's public healthcare system DATASUS^\ Here we 
consider that cities are the smallest administrative units with local government. The 
data consist of the number of homicides, suicides, and deaths in traffic accidents as 
well as the population for each city. 

US Data. We used data from the National Cancer Institute SEER, Surveillance 
Epidemiology and End Results downloaded from http://seer.cancer.gov/data/. The 
Institute provides mortality data aggregated for three or five years. To compare with 
Brazil, we use suicide mortality data for each American county and MSA 
(Metropolitan Statistical Area) accumulated through the five years of 2003 to 2007. 
Since the population is almost constant for a five years period, we adopted the average 
population as a measurement for the population on the allometry relation. 



Exponent Determination. In order to reduce fluctuations, we first apply a non- 
parametric fitting method to the data, namely, the Nadaraya-Watson kernel 
regression^^'^*^. We also compute the 95% (a = 0.05) confidence interval (CI) by the 
so-called a/2 quantile function over 300 random bootstrapping samples. The power- 
law exponent is calculated by performing an ordinary least square (OLS) fitting^^ over 
the results of the Nadaraya-Watson kernel regression in the population interval [10^, 
10*^] for Brazilian cities, and [10^ 10^ for American MSAs. 

Year 

Results Figure 2 | Temporal evolution of allometric exponent pior homicides 

Allometry in Urban Indicators. The main goal here is to investigate (red squares), deaths in traffic accidents (blue circles), and suicides (green 

the scaling behavior with city population of three urban indicators, diamonds). Time evolution of the power-law exponent P for each 

namely the number of homicides, deaths in traffic accidents and behavioral urban indicator in Brazil from 1992 to 2009. We can see that the 

suicides. For this, we analyzed data available for all Brazilian cities non-linear behavior for homicides and suicides are robust for this 19 years 

as well as suicide for US counties and MSAs, as presented previously period, and for the traffic accidents the exponent remain close to 1.0. 
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Figure 3 | Scaling relationship between suicides and population for US 
counties and MSAs. The small circles show the total number of suicides 
over five years (2003 to 2007) vs the average population for counties (a) 
and MSAs (b). The solid gray line indicate the best fit of a power law, using 
OLS regression, between the average total number of suicides and 
population. The dashed black lines delimit the 95% confidence band given 
by the Nadaraya-Watson kernel regression (solid black line)^^'^^. The 
allometric exponents are obtained through an ordinary least-squares 
(OLS) fit^*^ over the Nadaraya-Watson kernel regression applied to the 
suicides data. The values of the Pearson correlation coefficients p 
associated with these scaling relations are shown in each plot. We find ^ = 
0.87 ± 0.01 for counties and ^ = 0.88 ± 0.01 for MSAs with a 95% 
confidence interval estimated by bootstrap. The insets in each graph show 
the systematic decreases of suicide rates with population in both cases. 

on Materials and Methods. For 2009 in Brazil, as shown in Fig. 1, the 
increase in the number of casualties D with city population P for the 
three death causes can be properly described in terms of power-laws, 
D = DqP^, where Dq is a costant pre-factor, and the exponent P 
reflects a global property at play across the urban system. 
Interestingly, while the number of deaths by traffic accidents 
display isometric scaling, j5 ~ 1, homicides and suicides are both 
allometric, but obeying superlinear and sublinear scaling with 
population, respectively. These results suggest that the decision to 
commit a crime or to suicide, instead of being purely a consequence 
of individual choices, might have strong correlations with the 
underlying complex social organization and interactions. This does 
not seem to be the case of traffic accidents, since the strong evidence 
for isometric scaling, P = 0.99 ± 0.02, indicates that such events 
should result from random processes, i.e., no social relations need to 
be implied among the involved people. 

From the superlinear scaling exponent found for the number of 
homicides in Brazil at the year 2009, P = 1.24 ± 0.01, one should 
expect that, by doubling the population of a city, the number of 
homicides would grow approximately by a factor of 135% in average. 



instead of just growing 100%, as if we had an isometric scaling. This 
super-linear behavior is consistent with that found for serious crimes 
in USA^. The result obtained for suicide scaling seems quite surpris- 
ing. As depicted in Fig. Ic, the number of suicides scales with the 
population size as a power law with exponent = 0.84 ± 0.02, which 
implies an "economy of scale" of 22% in comparison with isometric 
scaling, similar to Kleiber's law for metabolic rates and animal 
masses. This sublinear behaviour is reminiscent of the seminal study 
by Emile Durheim\ one of the fathers of modern sociology. In his 
book Suicide, Durkheim explored the differences in average suicide 
rates among Protestants and Catholics, arguing that stronger social 
control among Catholics leads to lower suicide rates. The crucial 
contribution from Durkheim was certainly to treat the suicide as a 
social fact, by explaining variations in its rate at a macro level as a 
direct consequence of society-scale phenomena, such as lack of con- 
nections between people (group attachment) and lack of regulations 
of behavior, rather than individuals' feelings and motivations. 

The discrepancies observed in the scaling behaviors of homicides, 
deaths in traffic accidents and suicides become even more evident if 
we plot the average number of deaths per capita against city popu- 
lation, as shown in Fig. Id. Under this framework, the systematic 
decrease in suicide rate with population indicates that a large supply 
of potential social contacts and interactions might work as an "anti- 
dote" for this tragic event. This result is consistent with the idea that 
human happiness is more a collective phenomenon than a con- 
sequence of individual well-being conditions. In analogy with health, 
it is then possible to consider a "happiness epidemy" spreading in a 
social network^^. In Fig. 2 we show the dependence on time of the 
exponent p for a period of 18 years, from 1992 to 2009. We see a 
robust behavior for in such a way that , for different years, we still 
observe P > 1.0 for homicides, P < l.O for suicides, and P slightly 
above 1.0 for deaths by traffic accidents. 

We also analyzed data available for suicides in all US counties and 
MSAs. The data is an accumulation of the total number of suicides 
during a period of five years, from 2003 to 2007. In Figs. 3a and 3b we 
show the dependence of the total number of suicides during these five 
years on the average population for each county and MSA, respectively. 
As depicted, the number of suicides also scales with a sub -linear power 
law with exponent P = 0.87 ± 0.01 for counties, and P = 0.88 ± 0.01 
for MSAs, which are in agreement with our previous results for Brazil. 

In Fig. 4 we show the countour plot for the conditional probability 
of the rate of fatality, given a population size. In order to obtain the 
approximate density, we perform a kernel-density estimation over 
the sample of all Brazilian cities in 2009. The lines in Fig. 4 indicate 
the limits below which 10%, 50% and 90% of the data points are 
located. Besides confirming the superlinear, linear, and sublinear 
behaviors, these results also show how the probability distribution 
of rates of fatality varies with the population size. Also, the 10% and 
90% lines are representative of expected extreme cases of low and 
high fatalities, respectively. 

Discussion 

Allometric relations are ubiquitous in Nature, appearing in a wide 
variety of biological, sociological, chemical and physical systems^'^'^l 
Even the arrangement of Lego pieces has been recently reported to 
obey a sublinear scaling between the number of pieces types vs. the 
total number of pieces for many Lego sets^^. As originally proposed in 
biology by Kleiber^ and extended later by others^^, a single power law 
comprising 27 orders of magnitude can associate the metabolisms of 
a microscopical bacteria and a blue whale, weighting a few picograms 
and more than 100 t, respectively. Unlike allometric relations in 
Biology, few attempts have been made to explain the origin of uni- 
versal scaling laws describing urban indicators. In a recent study^^ a 
quantitative framework has been developed to consider the interac- 
tions between people over a social network that is capable to predict 
the allometric scaling for urban systems. 
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Figure 4 | Fatality per capita versus population for homicides, traffic accidents, and suicides. The color map represents the conditional probability 
density obtained by kernel density estimation. The bottom and top lines correspond to the 10% and 90% bounds of the distribution for each population 
size, that is 80% of the sampled points are between these lines. The middle line is the 50% level or "median" expected for each population size. The 
diagonal shape observed in the left side of density maps are cases of low number of fatal events, one or two fataUties. After this region we observe that the 
three level lines wiggle around an average power-law behavior. In the case of homicides the three level lines indicate an increase in the expected density of 
fatality with the population size. Similarly, for traffic accidents the lines are close to horizontal, that is, the probability distribution for the rate of fatality is 
near independent of the population. For suicides, the median shows a slight decrease with population size, while the 90% level, that is associated with cases 
of extreme rates of suicide, shows a pronounced decrease. The sublinear growth observed for suicides, as depicted in Fig. Ic, is likely due to the suppression 
of these extremely high rates in large urban areas. 



Here we studied the scaling relation with population of three 
behavioral urban indicators, namely, number of homicides, victims 
of vehicle crashes, and suicide events. We show that, unlike the 
incidence of vehicle crashes with fatal victims, which exhibits iso- 
metric (linear) scaling, homicides and suicides are characterized by 
allometric behaviors, superlinear and sublinear, respectively. 

Precisely, we found that the superlinear scaling relation for the 
number of homicides has an allometric exponent that lies between 
1.21 and 1.31 along the years, between 1992 and 2009. Despite all 
their positive attributes''^ (e.g., higher incomes and levels of creative 
activities), these results show that larger and urbanized cities also 
have a dark side, namely, higher levels of violence. Interestingly, the 
effect is exactly the opposite in the number of suicide events, which 
typically follow a sublinear scaling with population, with an allo- 
metric exponent value that varies in the range between 0.78 and 

0. 95. Such a result led us to conjecture that a suicide event should 
not be taken as an isolated individual decision. Many other hypo- 
thesis could be proposed to explain the observed higher rate of sui- 
cides of small cities. In any case, our idea is consistent with the 
conceptual framework put forward by Durkheim\ under which sui- 
cides need to be treated as social facts, actually affected by complex 
human relations. Our hypothesis is also in agreement with the recent 
finding by Schlapfer et al.^° that both the number of social contacts 
and the total communication intensity grow superlinearly with city 
population. The sublinear behavior found can be even considered as 
a counter- intuitive result. The more common and straightforward 
view would be to associate the suicide causes uniquely to a health 
condition of psychological or mental illness that could nevertheless 
be strongly linked to external urban factors. For instance, traffic jams, 
pollution, and stressful jobs all create a harmful environment for the 
population. However, the fact that we found sublinearity, namely, 
that the number of suicides is disproportionately small for larger 
cities, discloses an entirely different perspective. We conjecture that 
this property can be intimately related with an "emotional epidemy" 
as previously hypothesized^^'^^. This phenomenon can explain the 
systematic attenuation through the social network of contagious 
emotions states of potentially suicidal individuals. 
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